Formation and Interaction of Membrane Tubes 
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We show that the formation of membrane tubes (or membrane tethers), which is a crucial step in many 
biological processes, is highly non-trivial and involves first order shape transitions. The force exerted by an 
emerging tube is a non-monotonic function of its length. We point out that tubes attract each other, which 
eventually leads to their coalescence. We also show that detached tubes behave like semiflexible filaments with 
a rather short persistence length. We suggest that these properties play an important role in the formation and 
structure of tubular organelles. 
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PACS numbers: 87.16.Dg, 82.70.-y, 87.68,+z 

Biological membranes (such as the endoplasmic reticulum, 
the Golgi apparatus, the inner mitochondrial membrane, or the 
plasma membrane) often form highly dynamic tubular net- 
works [jl[]. The formation and transport of membrane tubes 
(with tens of nanometers in diameter) are thought to involve 
motor proteins that are able to grab the membrane and pull 
on it as they move along the filaments of the cytoskeleton 
[Oh. Nanotubes (or tethers) can also be pulled out by vari- 
ous experimental techniques (such as hydrodynamic flow [j^], 
micropipettes [Q], or optical tweezers [|]|), and very recently 
Roux et al. ^ have managed to set up a minimal in vitro ex- 
perimental system, in which tubes have been pulled by kinesin 
motor proteins. The physics of long tubular membranes is rel- 
atively simple and is well understood. However, the initial 
formation of nanotubes from planar or large spherical pieces 
of membrane is a subtle process, the understanding of which 
is crucial to the study of various biological processes (involv- 
ing tube formation or simply membrane pulling). In addition 
to tube formation we also study the interaction of membrane 
tubes pulled from the same membrane, and make a few com- 
ments on the mechanical properties of detached membrane 
tubes. 

Living cells maintain the surface tension of most of their 
membranes at a constant level by keeping lipid reservoirs (at 
fixed chemical potentials) p]. A constant pressure in closed 
vesicles is also maintained via osmosis. Therefore, for our 
study we choose the ensemble in which the surface tension 
a and the inside pressure p (relative to the outside) are fixed 
rather than the surface area of the membrane A or the volume 
of the vesicle V. 

With the bending term included the free energy of the mem- 
brane can be written as [f7|] 



T= J ~(2H) 2 dA + <jA-pV - fL 



(1) 



where k is the bending rigidity and H is the mean curvature 
of the membrane. The membrane is pulled in the Z direc- 
tion with a point force /, and the end-to-end distance of the 
membrane in this direction is denoted by L. 

For a tube of length L and radius R the free energy (at 



p = Q) can be written as ^tubc = [«/ (2i? 2 ) + a]2irRL — fL. 
To minimize J-'tube the surface tension acts to reduce the ra- 
dius, while the bending rigidity works against this. The bal- 
ance between the two sets the equilibrium radius R and force 
/o, which can be calculated by taking dFtube/dR = and 
dFtubc/dL = 0: 



Ro = 



(2) 



For typical values of k w 40 pNnm and a m 0.05 pN/nm one 
finds R m 20 nm and / w 12.6 pN. 

For simplicity we consider only axisymmetric surfaces with 
the Z axis being the symmetry axis (Fig. [j]a upper inset). Such 
surfaces can be parametrized by the angle ip(S), where S is 
the arclength along the contour. The coordinates R(S) and 
Z(S) depend on i/j(S) through 

R = cos'0 and Z = — smtp , (3) 

and the mean curvature can be expressed as 

2H = tp + (smtp)/R. (4) 

For axisymmetric surfaces one can derive the so-called gen- 
eral shape equation from the free energy ([[]) by variational 
methods [|[ \ 
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where a = a/ k = l/(2i?g) and p = p/n. Taking the first 
integral of this equation leads to [pfl 
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+cr sin ip - -pR - — , (6) 
with the integration constant / = f /(2ttk) = //(/o-Ro)- 
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For nearly flat membranes (ip <C 1), after the parameter 
change [ip(S), R(S)] — > ip(R), one can expand the shape 
equation (^) in powers of ip and keep only the terms up to 
linear order: 

R 2 ip" + Rip' - (R 2 d- + l)V = -fR - pR 3 /2 . (7) 

Because of the parameter change, the primes denote deriva- 
tions with respect to R (which, in linear order of ip, are iden- 
tical to derivations with respect to S). The general solution of 
this differential equation is 

V-OR) - i -- 5 + ^R + c 1 i 1 {rV&) + c 2 k 1 (r^) , (8) 

a R 2 a 

where L (x) and Ki (x) are modified Bessel functions, and ci 
and C2 are integration constants. l\(R\fa~') diverges expo- 
nentially for R — + oo. Because for a big vesicle we expect 
the shape to converge to that of a sphere, c\ must vanish. At 
R = the divergence of the 1 [R term must be canceled by the 
K\ term, leading to c 2 = —f/\fB. Integrating —ip(R) with 
respect to R gives the shape of the membrane in this linear 
approximation: 

(9) 

where the integration constant Zq serves as a reference coor- 
dinate, and we have expressed a, k, and p in terms of Ro, fo, 
and vesicle radius R vcs = 2a /p. The last term is a trivial con- 
tribution, describing a spherical vesicle under tension a and 
pressure p. The second term, which is proportional to /, is the 
linear response, and describes the deformation of the vesicle. 
The quantity between the brackets converges to [In (2) — 7] 
for R — > 0, where 7 = 0.577... is the Euler constant. For 
p = and large R the logarithmic term dominates, which cor- 
responds to a catenoid, the well known shape of a soap film in 
cylindrical geometry under zero pressure. 

Because the pressure makes only a trivial contribution, and 
has a negligible effect on tube formation from big vesicles 
(plug pR rj pRo — 2a Rq / R vcs <C a into the shape equa- 
tions), we neglect it from now on, and consider a piece of 
(initially flat) membrane that spans a ring of radius i? r i ng lo- 
cated at Z = 0. 

For large deformations (or pulling forces) the linear approx- 
imation breaks down, and we have to solve Eq. (|5]) together 
with Eq. (^J) numerically. Note that Eq. (||) could also be 
used, but it is numerically less stable. We start solving the 
differential equations from the ring, where we impose a zero 
curvature (or free hinge) boundary condition. Thus, the four 
initial parameters at Z = are as follows: (i) R = R T i ng ; (ii) 
ip = arcsin[//(27rcri? r i ng )] — e, where the small deviation e 
from the catenoid shape is chosen (with a shooting and match- 
ing technique) such that the contour line reaches the Z axis; 
(iii) ip — —(sinip)/R ensures zero mean curvature [see Eq. 
(fl)]; (iv) finally, ip is determined from Eq. (||). 

The results of the numerical solution can be seen in Fig. |IJ 
The main part of Fig. |I]a shows the shape (contour line) of the 
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FIG. 1: (a) The shape of an emerging tube at various lengths. The 
upper inset shows the parametrization of the surface, and the lower 
one illustrates the definition of L tu bc, which is the size of the de- 
viation of the shape from the linear approximation (dashed line), 
(b) Force vs. length curves for three different ring sizes. The inset 
shows the /-Ltube curve (solid line) and its asymptotic fit by Eq. 
^ (dashed line). 

membrane for i? r i ng = 20i?o an d different values of L. For 
small deformations (left line) Eq. (^|) gives a very good ap- 
proximation of the shape. In this linear regime the size of the 
deformation is approximately Ln n = Zi in (0) — Zn n ( i? r i ng ). 

For larger deformations a tube emerges in the middle, and 
the linear approximation fails. However, far from the tubular 
part, where ip is small, the approximation is still valid (Fig. |l]a 
lower inset). Thus, it is convenient to define the tube as the 
piece between Lu n and L, and the base as the rest of the mem- 
brane between and Lu n . This way, the dependence of the 
size of the total deformation L on the ring radius i? r i ng is ab- 
sorbed in the size of the base L\[ n , and the length of the tube 
^tube = L — L\i n (i.e., the deviation from the linear approxi- 
mation) becomes independent of i? r i ng . 

With this definition it is enough to determine the f-L curve 
for one particular ring size (e.g., 20i?o), from which the uni- 
versal (ring size independent) /--Ltube curve can be calcu- 
lated (Fig. |ljb inset). Because only the f-L curves have real 
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physical meaning, one can easily calculate them from the /- 
^tubc curve for any ring radius (3> Rq) by simply adding 
iiin(/, -Rring) to itube- The most intriguing feature of the /- 
L curves (shown in Fig. |ljb for three different values of i? r i ng ) 
is their non-monotonicity. The force first grows linearly, in 
accordance with the linear approximation (Q), and converges 
to fo for large L. But in between it overshoots by ~ 0.13/o, 
and then oscillates about fo with an exponentially decaying 
amplitude. This oscillation results in infinitely many intervals 
with negative slopes, which, in the /-ensemble (where / is the 
control parameter rather than L), are mechanically unstable 
and represent an infinite series of first order shape transitions 
at /o. Because Lu n is a monotonically increasing function 
of both / and i? r ing, the main peak becomes an overhang for 
large rings (> 20000i?o), leading to a first order transition 
even in the L-ensemble. 

For a nearly cylindrical section, like the tubular part of 
the membrane, the shape equations can be expanded in pow- 
ers of U(Z) = R(Z) — Rq. After changing parameters 
[tp{S), Z(S)] — > ip(Z), the expansion of Eq. (q|) up to first 
order in U(Z) and forp = reduces to 



Riu"" = -u 



(10) 



Here the primes denote derivations with respect to Z. The 
solution of this equation is the sum of two exponentially de- 
caying oscillations from the two ends (Lu n and L) of the tube: 
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where the integration constants converge to a\ ps 0.746, 
a 2 ~ 0.726, ct\ 0.347, and a 2 ~ 3.691 (determined by 
numerical fitting) as the length of the tube increases. These 
oscillations can be seen in Fig. || as overshootings from both 
the base and the tip of the tube. They can be understood in- 
tuitively by noticing that the smaller (larger) mean curvature 
near the base (tip) makes the effect of the bending rigidity on 
the tube radius less (more) pronounced, leading to a smaller 
(larger) R. 

The expansion of Eq. (^) up to the first non-vanishing order 
in U (Z) gives us the deviation of the force: 



f —A = Rtu'u'" - \rIu" 2 + l-Lu 2 



2Rl 



(12) 



Plugging Eq. (|TTJ) into Eq. ( |l2| ) yields 



/ — /o / -^tubc 

= 2ai<22 exp — ■ 



fo 



( ^tubc 

V2RoJ\y/2Ro 



Oil + a 2 



(13) 

which is an exponentially decaying oscillation as a function 
of Ltubo- This explains the observed non-monotonic behavior 
of the f-L curves, and gives a very good fit to the f-L tu b c 
curve (Fig. |ljb inset) with the same values of <Zj and on as 
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FIG. 2: The shape of a tube magnified in the radial direction (solid 
line), and the surface energy density (dashed line). 

determined above for the shape oscillations. In a different en- 
semble, Heinrich et al. [|llj also observed similar oscillations 
for axially strained vesicles. 

The surface energy density p = (k/2)(2H) 2 + a of the 
membrane is plotted in Fig. |[ It grows from a to 2a as we 
enter the tube from the base, and diverges at the tip. Because 
at the tip the angle ip is small again, Eq. (^|) can be used to 
describe the surface there, leading to a logarithmic divergence 
of the mean curvature [ |To[ 

In experiments, bundles of tubes can often be observed [^]. 
To study the interaction between two tubes, let us start with a 
planar membrane that spans a ring and is pulled perpendicu- 
larly by two point forces (fi and / 2 ) at a distance d. If d is 
large enough (R <C d i? r ; ng ), both protrusions (except 
for the vicinity of the points of pulling) can be described by 
the leading logarithmic term of the linear approximation (^). 
Their superposition can be used to calculate the <i-dependence 
of the free energy of the membrane. Note that fL in Eq. ([[]) 
has to be replaced by fiLi + f 2 L2, where L\ and L 2 are the 
sizes of the deformations. Straightforward calculation results 
in an attractive potential T(d) = const + 2Rq ln((i)/i/2//o 
between the two deformations. The same attraction has been 
found between membrane-bound adhesion molecules [JjjJ]. 
For the attraction between two tubes replace both /i and f 2 

by fo- 

To see what happens when two tubes get close, we per- 
formed numerical energy minimization with the SURFACE 
EVOLVER program rfTJJ]. We defined the initial topology as a 
piece of membrane spanning a ring of radius 12i?o at Z = 0, 
with two cylindrical protrusions in the positive Z direction. 
We applied zero curvature boundary condition at the ring and 
reflecting boundary condition at Z = 24i?o (Fig. ||). To make 
sure that the two tubes do not coalesce we placed a very nar- 
row (<C Ro) cylindrical obstacle between the tubes, perpen- 
dicularly to them, and at a distance h from the ring (not shown 
in Fig. [|). The obstacle simply pinned down the middle of the 
membrane at Z = h. We then measured the pinning force 
/ p (exerted by the obstacle on the membrane) as a function 
of h. Fig. ^] shows that / p is always negative (i.e., pushes the 
membrane), meaning that there is no energy barrier against 
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FIG. 3: The pinning force f p as a function of h. It goes from —fo/2 
(because if two tubes are far apart their attraction can be eliminated 
by a — fo/2 point force applied halfway) to — /o (because for large 
h the junction is pulled by two tubes upward but by only one tube 
downward). The solid line is a fit by an exponentially decaying os- 
cillation. 

coalescence, it occurs smoothly. This is consistent with the 
experiments of Evans et al. [Q]. 

Tubes can be detached via fission and then transported in- 
dividually in the cell as long prolate-like vesicles (keeping 
their original radius Ro). Because of volume and area con- 
straints their energy contains only the bending term: T = 
J (k/2)(2H) 2 dA. Although with such a high area-to-volume 
ratio the energetically most favorable shape is the stomato- 
cyte, these prelates are metastable [O]. They are also very 
flexible. If we consider a prolate as a rod and bend it with 
a curvature C (<C l/Ro), its energy, for symmetry reasons, 
increases quadratically with C. This energy increase can be 
written as (k/2)(XC 2 )2ttR L, where L is the length of the 
prolate, and the factor A = ^sin 2 0) — 1/2 indicates that 
during the integration of the energy around a cross-section 
(parametrized by <fi), only the out-of -plane component of the 
bending counts. Thus, the bending stiffness of a prolate-like 
vesicle is kttRq, and its persistence length uttRq/ {k^T) is in 
the order of a few hundred nanometers. 

Although these prelates are thicker than the microtubules, 
their persistence length is much closer to that of a DNA. This 
is because their wall is a two-dimensional fluid. Due to this 
fluidity they do not even resist twisting, which makes them 
ideal semiflexible filaments. 

We close this Letter by discussing some of the biological 
consequences of our theoretical results. The major biological 
relevance of the first force peak of the f-L curves is that in 
order to form a tube, the motor proteins must be able to pro- 
vide a force that is 13% larger than what is needed to pull a 
long tube. At a pulling force of / = /o the peak corresponds 
to a 2.1k « 21k^T high energy barrier, which is practically 



insurmountable for such a big and slow object as a growing 
tube. Thus, tube formation works on an all-or-nothing basis: 
motors can pull out tubes only if they are strong enough to 
overcome the major force peak. 

Although at the tip of the tubes the size of the lipids (~ 
0.5 nm) represents a natural cutoff length-scale for the diver- 
gence of the energy density, it goes up to tens of a. In terms of 
membrane rupture, the energy density acts as an effective sur- 
face tension, meaning that the most likely place for rupture to 
occur is the tip of the tubes. So if biological systems want to 
avoid rupture, they either have to protect the tips or distribute 
the pulling forces at larger areas. That could be achieved, e.g., 
by utilizing cap proteins or lipid rafts. 

We have shown that without external pinning tubes coalesce 
smoothly. Thus, to explain bundle formation, other physical 
effects that could prevent the coalescence of tubes (e.g. adhe- 
sion between the tubes and the cytoskeleton) must be taken 
into account. 
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EU Marie Curie Fellowship HPMF-CT- 1999-00 124. 
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